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Lyman alpha emission from the first galaxies: Signatures of 
accretion and infall in the presence of line trapping 
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ABSTRACT 

The formation of the first galaxies is accompanied by large accretion flows and viri- 
alization shocks, during which the gas is shock-heated to temperatures of ~ 10^ K, 
leading to potentially strong fluxes in the Lyman alpha line. Indeed, a number of 
Lyman alpha blobs has been detected at high redshift. In this letter, we explore the 
origin of such Lyman alpha emission using cosmological hydrodynamical simulations 
that include a detailed model of atomic hydrogen as a multi-level atom and the effects 
of line trapping with the adaptive mesh refinement code FLASH. We see that baryons 
fall into the center of a halo through cold streams of gas, giving rise to a Lyman alpha 
luminosity of at least 10^^ erg s~^ at z = 4.7, similar to observed Lyman alpha blobs. 
We find that a Lyman alpha flux of 5.0 x 10~^^erg cm~^ s~^ emerges from the enve- 
lope of the halo rather than its center, where the photons are efficiently trapped. Such 
emission can be probed in detail with the upcoming James Webb Space Telescope 
(JWST) and will constitute an important probe of gas infall and accretion. 
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1 INTRODUCTION 

Early- type galaxies are known to produce copious Lyman al- 
pha photons due to their spatially extended gas distribution. 
In protogalactic halos gas transfers its gravitational binding 
energy to the excitation of hydrogen atoms which results in 
Lyman alpha emission (|Haiman et al ] l200Ql :[ Diikstra et al.l 
l2006al lbh. Therefore, the gas in dark matter halos exceed- 
ing a virial temperature of ~ 10^ K may be detected in 
the Lyman alpha line emission. The presence of ionizing ra- 
diation sources may enhance the Lyman alpha flux as gas 
is photo- ionized, and and supernova feedback may further 
increase the escape fraction by generating a more clumpy 
structure. These ionization sources could be the first stars 
or mini-quasars. A luminous quasar can boost the emis- 
sion of Lyman alpha ph otons by several orders of magnitude 
iHaiman k Ree3 (l200lh . 



Many Lyman alp ha blobs (LABs) have been ob- 
served at high redshi f t (I Yang et al" l2QQ9l: jouchi et al.ll20Q9l : 
iMatsuda et aD l2004l : ISteidel et al.1 |20QQ|). In the light of 
recent detections a t redshift > 7 (jLehnert et al.l 2010l : 
IVanzella et al.l 1 20 id ) , it is of high interest to understand 
what drives the emission and how it is spatially distributed. 
A number of LABs have been observed whose most prob- 



able origin is cold accretion of gas onto d ark matter halos 
(iNilsson et al.ll2006l : ISmith Jarvisll2007l ). 

Numerical simulations show that cooling by Lyman 
alpha radiation induces collapse in protogalactic halos. 
Moreover, baryons accumulate into the center of ha- 
los by penetration of cold streams of gas through the 
shock heated medium jPardal et al.l 200ll : Wise et al. 20081: 
IJappsen et al.' '200 91: iRegan k Haehneltl l2009l: IPekel et al.l 
2009; Keres et al. 2 Q09| : lGoerdt et al.ll2010l : lLatif et al.ll2010l : 
[Shang et al.ll2010l ). Cold streams with temperatures of ^ 
10^ K could be pote ntial sources of spatial l y extended Ly- 
man a lpha emission (jPiikstra k Loebll2009l l . Ijohnson et al.l 
(|2010l ) found that Lyman alpha radiation can also be emit- 
ted during accretion of gas on black holes formed by direct 
collapse in the first galaxies. Such emission is potentially de- 
tectable with JWSTIj. Similarly, it allows to probe the star- 
burst component th rough the enhanced emission in several 
recombination lines ([Johnson et al.|[2009l V 

The presence of large columns of neutra l hydrogen gas 
causes the trapping of Lyman alpha photons ([Spaans Silkl 
l2006l : iLatif et al.1 l2010l : ISchleicher et"all l2010l ). which was 
neglected in some of the previous studies. In this letter, 



http: / /www. stsci.edu/jwst/instruments/nirspec/sensitivity 
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we aim at a self-consistent modeling of Lyman alpha emis- 
sion driven by accretion flows, including dynamics, non- 
equilibrium chemistry of H, H^, He, He^ and He^^, de- 
tailed level populations of atomic hydrogen, as well as the 
trapping of hydrogen line photons due to large column den- 
sities. The prime objective of this work is to study the origin 
and spatial distribution of Lyman alpha emission, which is 
detectable with JWST. 



2 MODELING OF THE PHYSICS 

For our simulations, we employ an extended version of 
the adap tive mesh refinemen t (AMR) hydrodynamics code 
FLASH (jPubev et alJ l2009l ). FLASH is a module based 
Eulerian grid parallel simulations code. We use an AMR 
grid to achieve high dynamic resolution in the regions 
of interest. We employ the directionally split piece wise 
parabolic method (PPM) for hydrodynamic calculations, 
which is an improved fo rm of the Godunov method 
(jColella Woodward! Il984h . This method is well suited 
for flows involving shocks and contact discontinuities. The 
dark matter is simulated based on the particle mesh (PM) 
method. In order to create Gaussian random field initial 
conditions, we run the COSMICS package developed by 
lBertschinge3 (Il995l ). We start our simulation at redshift 
100. Our computational periodic box has a comoving size 
of 10 Mpc. We perform our simulations in accordance to 
the ACDM model with WMAP 5-years parameters. We se- 
lect a 4 X 10^ M0 halo at redshift 7 and follow its collapse. 
We enforce 8 additional levels of refinement, which gives 15 
levels of refinement in total. In this way, we obtain a dy- 
namic resolution limit of 70pc (co- moving), and can still 
follow the further evolution until redshift z = 4.7. We im- 
pose th e Truelove criterion b y refining according to the Jeans 
length (jTruelove et al.lll997l ). We resolve the Jeans length by 
at least 20 cells. Such resolution was shown to be sufficient 
for accurately resolving turbulent structures during gravita- 
tional collapse (Federrath et al. 2010). When the highest re- 
finement level is reached, we prevent artificial fragmentation 
via Jeans-heating. We assume a primordial gas composition 
with 75% hydrogen and 25% helium by mass. 

We have developed a chemical network for the non- 
equilibrium modeling of non-molecular species. For this pur- 
pose, we solve the rate equations of the following species 
H,H+,He,He+,He++ and e" for non-equilibrium ioniza- 
tion. The rate equations for these species are solved us- 
ing the backward d ifference formula (BDF) method of 

rates of 
We fur- 



Anninos et al. (| 19971 ) . We adopt the chemica l 



Abel et al .1 (|l997f ) and ISchleicher et"all dioQ; 



ther solve for the level populations of atomic hydrogen 
up to the fifth electronically excited state, and model 
the non-equilibrium cooling including hydrogen line emis- 
sion, collisional ionization cooling, recombination cooling, 
Bremsstrahlung cooling and Compton cooling/heating. The 
tr ansition rates fo r the level populations are based on work 
of lOmukail (|200lh . At an optical depth tq > lO'^, the pho- 
ton escape time becomes longer than the gas free fall time. 
Because of the weak dependence of the photon escape time 
on the gas number density, tph oc n~^/^, trapping becomes 
important important during the collapse since the latter 
scales as n~^/^. A detailed discussion of such line trapping 



effects is given by lOmukail (l200lh : ISuaans k Sil^ (l2006l ): 
ISchleicher et"all (l2010l ). We have computed the Lyman al- 
pha luminosity by determining the Lyman alpha emissivity 
and integrating it over the virial volume of a halo. Further 
de tails of the luminosity and flux calculation can be found 
in lDiikstra et al.1 (l2006a 



3 RESULTS AND CONCLUSIONS 

The initial density perturbations decouple from the Hubble 
flow and begin to collapse via gravitational instabilities. The 
gas is shock-heated during the non-linear collapse of density 
perturbations. At redshift 10, the halo begins to virialize. In 
the process of virialization, part of its gravitational energy 
is converted into thermal energy. This heats up the gas and 
results in the emission of Lyman alpha photons. The gas is 
accreted into the center of the halo through cold streams 
of ^ lO^K as shown in figure [T] The upper panels of figure 
[T]show the density at different redshifts while lower panels 
show corresponding temperatures. Typical number densities 
of the cold streams are of the order of 0.01-1 cm~^. Density 
and temperature radial profiles are depicted in the upper 
left and right panels of figure [21 We see that in the presence 
of hydrogen line emission the gas cools isothermally. The 
density profile of the halo is ^ r~^"^, which agrees with re- 
sults of I Wise et aP (|2008l ). The ionization degree of the gas 
is shown in the bottom left panel (HII abundance) of figure 
[21 We see that most of the gas remains neutral at 10^ K, ion- 
ized fraction goes down due to faster recombination at higher 
densities ^ n~°"^. The gas radial velocity is depicted in the 
bottom right panel of figure [21 It can be seen from the figure 
that gas is falling into the center of the halo. The velocity 
profile also shows that gas is accreted onto the halo through 
accretion shocks. The cold streams have typical column den- 
sities of about lO^^cm"^. A radial profile of the column den- 
sity is shown in figure [3l At columns above lO^^cm"^ the 
gas optical depth exceeds 10^ and Lyman alpha trapping 
becomes effective. The Lyman alpha emissivity is shown in 
figure [3l It can be seen from the figure that above columns 
of the order of lO^^cm"^, Lyman alpha photons are trapped 
and cooling through this line is completely suppressed. That 
is why the Lyman alpha emissivity sharply declines in fig- 
ure [3l We explored the local density variations and found 
that the maximum variation of the column density does not 
exceed an order of magnitude. In this sense, our calculation 
should provide a conservative lower limit on the expected 
flux. In the presence of Lyman alpha trapping, cooling can 
proceed through highe r electronic states of atomic hydrogen 
(ISchleicher et al.ll2010l ). For comparison, the total emissiv- 
ity from higher states of atomic hydrogen and recombina- 
tion/Bremsstrahlung processes is shown in figure [3l The to- 
tal emissivity plot shows that despite line trapping, cooling 
still proceeds through these higher states of atomic hydro- 
gen, particularly through 2s- Is and 3-2 transitions, and the 
thermal evolution is approximately isothermal. 

The radial profile of the enclosed Lyman alpha luminos- 
ity is shown in figure [3l If the emission originates preferen- 
tially from the center, an approximately flat profile would be 
expected, while for uniform emission, a power-law behavior 
as ~ would be expected due to volume effects. Here, we 
find the absence of emission in the central core, then a sharp 
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Figure 1. The upper panels of this figure show the density shces through the center of the halo. The left panel is for redshift 4.7 and right 
panel for redshift 8. Velocity vectors are overplotted on the density slices. The bottom two panels show the corresponding temperatures. 
The figure shows the inner region of 20 kpc in comoving units. 



by four orders of magnitude, 



and 



more modest increase 
between radii of 10"^^ cm — 10^^ cm. This behavior reflects 
the generation of Lyman alpha emission through accretion 
shocks. JWST can confirm the presence of such a brightness 
profile. The total Lyman alpha luminosity from the halo is 
10^^ er g s~^, which is consiste nt with ob served Lyman alpha 
blobs (lOuchi et a l. 2009; Goer dt et all EoiO). 

The emerging flux from the halo is depicted in figure 
[3 We find a total flux of 5.0 x 10"^'^ erg cm"^ s"^ At red- 
shift 4.7, the observed Lyman alpha wavelength is at 0.68 
microns. This can be detected with the JWST-instrument 
NIRSpec for integration times of 10^ s, with S/N=10 and 
R=100. JWST spectrograph NIRSpec will have angular res- 
olution of 0.1 arc-sec for 2 micron, and NIRcam will be well- 
suited for higher angular resolution with field of view of 
2.2x2.2 arcmin^ and 4096^ pixels for shorter wavelength of 
0.6-2.3 micro meter. The halo in our case will have an angu- 
lar size of 0.5 arc-sec at redshift 4.7. It should be detectable 
with NIRcam/NIRSpec. For extended emission of Lyman 
alpha, it may need higher integration time 10^ sec) to 
resolve the flux from extended sources. The total mass of 
our halo at redshift 4.7 is 2 x lO^^M©, which is not extreme 
in any way. Selecting higher mass halos will produce higher 
fluxes as there is a power law relation between the mass of 



a halo and its lumin osity Lly^ oc M^^^ (jPekel et al.l l2009l : 
iDiikstra et "alll2006al ). We stopped our simulation at a red- 
shift of 4.7 as it becomes computationally too expensive to 
follow the further time evolution. 

In this work, we assume that there is no X-ray/UV 
background flux. Its presence c an heat the ^as and may 
increase the Lyman alpha flux f Spaans Meiierinkl l2008l : 
IZaroubi et al. 12007 ). We have compared our results with 
(IDiikstra et al.ll2006al ) and found good agreement. We have 
assumed here that the halo is metal free. The addition of 
small amounts of dust can absorb Lyman alpha photons effi- 
cient ly, unless the gas is inhomogeneous (|Haiman Spaana 
Il999l ). However, in- falling gas may still be pristine for iso- 
lated halos. We have ignored H2 cooling, which is suppressed 
in the pre sence of a strong UV background (J21 > 100) 



Diikstra e t al. 2004; Cazaux & Spaans 2 0091 ) . Ishang et al.1 
2OIOI ) found that even J21 30 will be suffi cient to quench 



H2 for mation for local variations in flux see IDiikstra et al.l 
(|2008l V At higher redshifts, it is possible that H2 may form 
for more modest radiation backgrounds. Even then, we ex- 
pect H2 formation preferentially at the center of the halo, 
while the accretion at the virial radius may still include a 
gas phase at 10^ K. This supports our main conclusion that 
Lyman alpha radiation originates mostly from gas in halo 
envelope. 
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Figure 2. All panels in this figure are log log plots. The upper left panel shows the radial density profile of the halo. The radial 
temperature profile for the halo is depicted in the upper right panel. HII abundance is shown in the lower left panel. The radial velocity 
of the halo is depicted in the lower right panel. 
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Figure 3. All panels in this figure are log log plots. The upper left panel shows column density plotted against the radius of the halo. 
Lyman alpha emissivity radial profile emerging from the halo is shown in the upper right panel. Lyman alpha luminosity is shown in the 
lower left panel. Lyman alpha flux is plotted against the radius of the halo in the lower right panel. 
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Some L ABs are also associate d with massive star form- 
ing galaxies (|Matsuda et al.|[2006l l where stellar feedback or 
an AGN could power Lyman alpha radiation. It is clear that 
starbursts and AGNs are also potential candidates to power 
or enhance the observed Lyman alpha luminosities. It is not 
fully clear if our results are also applicable to such situa- 
tions, although we expect similar effects of line trapping in 
the centers of such galaxies if Nh/Av ^ 10^^cm~^km~^s, 
with Av the cumulative velocity difference along the atomic 
hydrogen column Nr. Current modeling uncertainties con- 
cern both the star formation efficiency and the initial mass 
function. Both radiative feedback, leading to the formation 
of HII regions, as well as mechanical feedback, providing a 
more clumpy structure with higher escape fractions, need 
to be modeled self-consistently. Outflows if present may en- 
hance the detectability of Lyman alpha emission in the out- 
skirt of the halo, 5 % of the emitted Lyman alpha photons 
could be directly transmitted to t he observer along the line 
of sight ([Piikstra Wvitli3l2010h . In the future, cosmolog- 
ical radiative transfer simulations including feedback effects 
should be performed to obtain more robust results. 
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